Interacting Quantum Dot Coupled to a Kondo Spin: A Universal Hamiltonian Study 
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We study a mesoscopic interacting quantum dot described by the "universal Hamiltonian" that 
is coupled to a Kondo spin. The ferromagnetic exchange interaction within the dot leads to a 
stepwise increase of the ground state spin; this Stoner staircase is modified non-trivially by the 
Kondo interaction. We find that the spin-transition steps move to lower values of the exchange 
coupling for weak Kondo interaction, but shift back up for sufficiently strong Kondo coupling. The 
problem is solved numerically by diagonalizing the system Hamiltonian in a customized good-spin 
basis and analytically in the weak and strong Kondo coupling limits. The interplay of Kondo and 
ferromagnetic exchange can be probed with experimentally tunable parameters. 
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A singly-occupied localized electron level (an "impu- 
rity" spin) interacting with a delocalized electron gas is a 
paradigmatic system in quantum many-body physics. It 
gives rise to the non-perturbative Kondo effect in which 
the localized electron's magnetic moment is fully screened 
by the delocalized electrons below the Kondo tempera- 
ture T/f. The Kondo effect was well-studied in the con- 
text of an impurity moment embedded in a bulk metal 
Recently, Kondo physics has been the subject of 
much renewed interest [2|, following its observation in 
experimentally tunable quantum dots [3|. 

These experimental advances have been accompanied 
by progress in the theoretical treatment of the mesoscopic 
Kondo problem [3, d, @, 0, H, 0, [IB] . In these mesoscopic 
systems, either the electrons in the leads or the electrons 
in a large dot play the role of the "electron gas" and a 
small spin- 1/2 dot represents the Kondo spin. Here we 
focus on the latter case [see Fig.QJa)], where the discrete 
single-particle level spacing 5 and mesoscopic fluctuations 
of the large dot may alter the Kondo effect when Tjc ~ 5 

UMMMM- 

A formidable challenge for both mesoscopic and bulk 
Kondo theory is to take into account electron-electron 
interactions in the electron gas. In the mesoscopic case 
this task is simplified when the electron gas is confined 
to a large quantum dot in which the electron dynamics 



is chaotic [llj, [12j]. When the dot's Thouless energy Et 
is large compared with S, the effects of electron-electron 
interaction are captured by the so-called universal Hamil- 
tonian (UH) [l3| . valid in an interval Et around the dot's 
Fermi energy Ep. For a fixed number of electrons, the 
dominant interaction in the UH is a ferromagnetic ex- 
change interaction. Detailed comparison between theory 
and experiment for the statistics of Coulomb blockade 
peak heights and spacings shows that including the UH 
ferromagnetic exchange term is both necessary and suf- 
ficient to obtain quantitative agreement [3]. Thus, one 
can use this UH to obtain an experimentally-relevant de- 



scription of a large interacting dot (henceforth called the 
"dot") that is Kondo-coupled to a small dot with odd 
electron occupancy (henceforth called the "Kondo spin" 
Sk)- Such a model was first discussed in the framework 
of a mean-field approximation 0], where Kondo corre- 
lations in a dot close to its ferromagnetic Stoner insta- 
bility, J s ~ S, were investigated. A regime just below 
the instability was identified where the Kondo coupling 
substantially reduces the dot's polarization. In contrast, 
studies in the bulk [3] found that a Kondo impurity en- 
hances the polarization of a surrounding gas of electrons 
at similar high values of J s < 5. 

The simplicity of the UH allows one to look for sig- 
natures of the Kondo interaction in the magnetic prop- 
erties of the dot without necessarily assuming that the 
dot is very close to its bulk Stoner transition (in a quan- 
tum dot, J s is typically a fraction of S 1J, LLJ]). Stan- 
dard numerical methods for the Kondo problem such as 
quantum Monte Carlo [t| and numerical renormalization 
group (NRG) techniques [10( are, however, not easily ap- 




FIG. 1: (Color online) (a) Schematic diagram of a small quan- 
tum dot with spin Sk (Kondo spin) that is coupled antiferro- 
magnetically (coupling constant Jj,) to a large quantum dot 
with spin S^. The large dot, described by the universal Hamil- 
tonian, is characterized by a ferromagnetic exchange interac- 
tion (coupling constant J s ). We assume the large dot to have 
N equally-spaced (spacing 8) single-particle levels in a band 
of width 2D (half- filling), (b) The large dot is represented in 
the site basis (squares), in which Sk couples only to site 0. 
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plied to this problem because the ferromagnetic exchange 
coupling introduces a sign problem in the former and 
non-local correlations in the latter. Here we use a cus- 
tomized diagonalization method that takes advantage of 
the global spin rotation invariance, using the good-spin 
eigenstates of the UH as a basis [lj, EH . 

The Hamiltonian of our system, schematically illus- 
trated in Fig. [TJ is given by d, EH , 

H= ]T e^cl a c^-J s S 2 d + J k S K - Sd (0). (1) 

/i,er— ± 

The first two terms in (p} constitute the UH of the dot 
[IH (ignoring a constant charging energy and a Cooper 
channel term), described by N spin-degenerate single- 
particle levels e M and spin S d = \ X^o-o-' c \ ia T <y<y' c ^<y 
(r are Pauli matrices). The coupling of the dot to the 
Kondo spin Sk (with Sk = 1/2) is mediated by its 
spin density s d (0) = ^E^' <M ) c j«x' r <T<T'</>* (0)cw at 
the tunneling position r = [^(r) is the orbital wave 
function of level p] . The parameters J s and J k are the 
exchange and Kondo coupling constants, respectively. 
In this study we ignore mesoscopic fluctuations, taking 
equally spaced single-particle levels [covering a band of 
width 2D = (N-l) x 8} and </y(0) = l/\fN. We also 
assume half-filling of the band so that the number of dot 
electrons is N. The average local density of states of the 
dot is p = 1/(NS) 0. 

The spin-rotation invariance of the Hamiltonian ([T]) 
implies the conservation of the total spin Stot = S^ + Sk , 
so that Stot and Stot = -^tot are good quantum num- 
bers. To take advantage of this symmetry, we construct 
a good total spin basis by coupling the eigenstates of 
the UH with those of the Kondo spin. The UH eigen- 
states with dot spin Sd are characterized by \~fSdMd) (7 
denotes orbital occupations and other quantum num- 
bers distinguishing between states of the same dot spin 
Sd)- Thus a basis of the coupled system with good total 
spin is \~f SdS tot M tot) (for simplicity the quantum number 
Sk = 1/2 is omitted). In this basis the UH is diagonal 
with energies ^2 e M n M — J s S d (S d + 1). The Kondo term 
Hk = Jk Sk • Srf(O) is a scalar product of vector opera- 
tors in the uncoupled spaces. Thus, its matrix elements 
in the coupled basis conserve Stot ; -^tot and are given by 

(yS' d S tot M tot \H K \iS d S tot M tot ) = J k (-l) s * +1 / 2 + s ^ 

in terms of a Wigner-6j symbol and the reduced matrix 
element of the spin density s d (0) (known in closed form 
[lit). ■"■ n this formulation, the full Hamiltonian H has a 
block diagonal structure in Stot,Mtot- 

The problem is further simplified by transforming to 
the basis of sites i (0 <i < N — 1), in which the one- 
body part of the dot's Hamiltonian is tridiagonal and 



Hk = Jk^K ■ so with s, the spin at site i [l[ (see Fig.QJi). 
The exchange interaction is invariant under such trans- 
formation and has the same form as in Eq. ([1]) with 
Sd = J^i^) 1 s i- We can thus recast our formalism in this 
site basis, where only neighboring sites are coupled and 
the Kondo spin interacts solely with site 0. Due to these 
features the many-body Hamiltonian matrix in the site 
basis is more sparse than in the orbital basis, allowing 
for an efficient diagonalization in each subspace of good 
Stot using a Lanczos-Arnoldi algorithm. In this approach 
we can conveniently diagonalize |T]) for dots with up to 
N ~ 12 levels, where the total Hilbert space contains 
~ 5.4 x 10 6 basis states. 

We calculated the lowest many-body energy eigenvalue 
for each value of Stot and thereby determined the ground- 
state value of the total spin for different values of J s 
and Jk- This quantity has been studied theoretically 
0, [i| and can be probed experimentally [13] . As J s in- 
creases, the ground state spin S to t is expected to un- 
dergo successive transitions to higher values (known as 
the Stoner staircase) until the dot becomes fully polar- 
ized at J s ~ 5. For Jk — » 0, the spin transitions oc- 
cur at J™ = 5(m + l)/(m + 2), with m = 1,3,5,... 
(to = 2, 4, 6, ... ) for an odd (even) number of dot elec- 
trons N. These transition steps in the Stoner stair- 
case are shifted by the Kondo interaction. In Fig. [2ji,b 
we show the transition curves (colored lines) separat- 
ing regions of fixed ground-state spin Stot in the two- 
dimensional parameter space of J s , Jk- We observe that 
these curves are monotonically decreasing for Jkp < 1 
and monotonically increasing for J^p > 1. In the strong- 
coupling limit, they converge to values of J s that are 
either lower (for smaller S to t) or higher (for larger S to t) 
than their corresponding weak-coupling values J™ . 

To gain insight into the behavior of the transition 
curves, we evaluate them for weak Kondo coupling in 
first-order (degenerate) perturbation theory. We find 

JT = [(to + 1)6 - a m J k ]/[m + 2] , (3) 

where a m are positive constants of order one. These 
perturbative results (dotted lines in Fig. [2ji,b) agree well 
with the numerical calculations for Jkp < 0.1. 

The negative slope of the transition curves at weak 
coupling can be understood by considering that in this 
weak-coupling regime the Kondo spin Sk plays the role 
of an effective magnetic field, polarizing the dot in the 
direction opposite to its own spin. This will favor larger 
values of S d and thus also larger values of S to t, hence the 
negative slope. As we approach the Stoner instability 
J s — * <5, (i.e., increasing to), the gain in exchange corre- 
lations J s S d (S d + 1) dominates over the gain in Kondo 
correlations, hence the flattening of the slope in this limit. 

A perturbative analysis can also be carried out for the 
limit of strong Kondo coupling, which, at zero temper- 
ature, is characterized by Tk 3> S. For dots with suf- 
ficiently large N, this limit can be reached already for 
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FIG. 2: (Color online) Ground-state spin Stot of the system in 
Fig. Q] at finite exchange J a and Kondo coupling for an odd 
(left column) and even (right column) number of electrons N. 
(a),(b) Transition curves for N = 11 (left) and AT = 12 (right), 
separating regions of fixed Stot- Numerical results (colored 
solid lines) are compared with analytical estimates in the weak 
and strong coupling limits (dotted lines). (c),(d) Spin transi- 
tion curves for fixed band width 2D but different N (top to 
bottom): N = 5,7,9, 11 (left) and N = 6, 8, 10, 12 (right). For 
increasing N (at JkP^$>l), the curves converge to the Stoner 
staircase of a dot with N — l electrons (colored arrows). 



J up -C 1, where the Kondo temperature is of the or- 
der T K ~ De~ l ' Jk P (at J s = 0) [1]. However, for the 
present case of N < 12, the limit Tr- ^> 5 requires 
JkP ^ 1 or, equivalently, Jt > 2D. The latter condi- 
tion represents the bare strong-coupling limit for which a 
perturbative solution is available (without renormalizing 
the band width) [l8j]. In this limit the Kondo spin Sk 
and the spin So at site form a strongly bound singlet 
Sk.o = {Sk,o = S_r- + s ) that is effectively decoupled 
from the rest of the spin-chain with sites i > 1. The 
tridiagonal one-body Hamiltonian of sites 1 < i < N — 1 
can be rediagonalized to give X^o- ^cj^c^, describing a 
"reduced" dot with new orbital wave functions (p^ (r) and 
single-particle energies e M . This dot has one less level and 
one less electron than the original dot. While the origi- 
nal dot levels are equally spaced, the level spacings in the 
reduced dot are given by e p+ i — = 5^ 8 + > 5, 

(/3 M > are of order 1 and increase monotonically from 
the new band center towards the band edges). 

To explore how the strong-coupling limit is modified 
in the presence of exchange interaction in the dot, we 
rewrite the latter as — J S S^ = — J S S^ — 2 J s s • — 2>J S /A 



where Sd = Sta 1 s « i s the spin of the reduced dot. The 
cross term — 2 J s s$ ■ Sd has vanishing matrix elements in 
the singlet subspace Sk.o = but induces virtual tran- 
sitions to the triplet subspace Sk,o = 1 that renormalizc 
the exchange coupling constant J s — ► J s = J s (l + J s /Jk) 
(details will be presented elsewhere). To lowest order in 
1/Jfc, our system is thus described by an effective Hamil- 
tonian ^4ffV _ ^ s ^ *hat ^ as ^ e f° rm 01 a UH for 
the reduced dot with single-particle energies and ex- 
change constant J s (J s — > J s for J^— >oo). The spin tran- 
sition curves of this reduced dot (dotted lines in Fig.[5h,b) 
are found to be in good agreement with the exact numer- 
ical curves when Jip> 1. Despite their accuracy in the 
appropriate limits, the perturbative estimates do fail in 
the fully Kondo correlated intermediate regime and our 
exact numerical solutions are essential to show that there 
is a smooth crossover between the limits. 

The spin transition curves in the crossover from weak 
to strong coupling (see Fig. ar e determined by two 
counteracting effects: (i) The effective removal of an elec- 
tron from the dot shifts down the Stoner staircase accord- 
ing to J™ — * J 7 / 1 ^ 1 - Since the reduced dot has one less 
electron, the shifted Stoner staircase is associated with 
the opposite number parity of electrons, (ii) The effective 
removal of a level from the dot stretches the step size in 
the staircase due to the larger level spacing in the reduced 
dot (i.e., S/i > S), and thus increases the spin-transition 
values of J s . The downward shift in (i) is independent 
of N, but weakens for increasing m (where the step val- 
ues J] 71 are more densely spaced). The upward shift in 
(ii) is a finite-size correction ~ 1/N that decreases with 
N, but increases with m because of the non-uniform 5^. 
For smaller values of J Sl effect (i) dominates over (ii), 
resulting in an overall downward shift of the transition 
values in the strong-coupling limit (as compared to the 
weak-coupling values J™). Close to the Stoner insta- 
bility, however, finite-size effects (ii) dominate over (i), 
leading to transition values larger than J]". 

To investigate the interplay between effects (i) and (ii) 
more closely, we compare the spin transition curves of our 
original systems (N = 11, 12) with systems of equal band 
width 2D, but different values of N. Results shown in 
Fig- Hfcjd demonstrate that finite-size effects (ii) decrease 
with increasing N, leading to a convergence of the strong 
coupling transition curves towards J™ -1 . 

Dots with a large band width (D 3> 5) could, in princi- 
ple, be studied by extending the NRG method to include 
an exchange interaction. Truncation of the band width 
below Tk leads to a strong-coupling effective Hamilto- 
nian that includes additional interaction terms. It would 
be interesting to investigate the effects of these terms. 

Signatures of the interplay between the intra-dot ex- 
change and the Kondo coupling are revealed by applying 
an in-plane field B [T3|, adding a Zeeman term gusB St ot 
to the Hamiltonian H in Eq. |T]) (g is the gyromagnetic 
factor and /ib the Bohr magneton). This term commutes 
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with H (although it breaks the M to t degeneracy) and fa- 
vors a parallel configuration of dot spin (|) and Kondo 
spin (fr), increasing S to t at 4 = 0. The addition of 
the Kondo interaction opposes such a parallel alignment 
(Tit); correspondingly we find (Fig. [3k.) the spin tran- 
sition values of B monotonically increasing vs. Jk- A 
more complex and subtle behavior appears for J s 
(see Fig.[3ti), where non-monotonic spin-transition curves 
arise. The behavior in weak Kondo coupling can again 
be understood in perturbation theory, for which the dot 
spin can still be regarded as a good quantum number. 
In the limit Jk — ► the ground state will always be 
the parallel configuration (TIT) aligned with the exter- 
nal field, so that the spin transition lines slowly increase 
with Jfe (as for J s = 0). At larger values of Jk the en- 
ergy of the anti-parallel (J, -ft-) configuration (with rear- 
ranged orbital occupancies) becomes lower in each spin 
subspace. This happens at lower Jk for the lower value 
of Stot, leading first to a marked increase in the slope 
of the transition line. Increasing Jk further one reaches 
the point at which the anti-parallel configuration is also 
favored in the subspace with higher Stot- This decreases 
the slope of the spin transition line, making it negative in 
some cases. For even larger Jk the perturbative picture 
breaks down and the transition curves make a smooth 
crossover to the strong coupling picture where the ef- 
fective exchange interaction constant J s now decreases 
with increasing Jk, favoring lower Stot again and giving 
a positive slope to the transition curves. Initial calcu- 
lations with non-uniformly spaced single-particle levels 
show that such a non-monotonic behavior persists in in- 
dividual dots. We defer detailed studies of mesoscopic 
fluctuations to future work. 

The ground state spin of quantum dots (in the absence 
of Kondo coupling) has been measured in a number of 
experiments [17|; typically this requires varying the in- 
plane (Zeeman) field B, allowing us to map out the spin 
transition diagram in the following manner. The points 
of degeneracy between ground states with different spin 
for fixed Jk are determinable by kinks in the slope of 
the Coulomb blockade peak positions vs. in-plane field B 
[lij ]. Tuning Jk (by means of a pinch-off gate) at fixed 
J s will cause the kinks to shift to higher or lower values 
of B in a manner predictable from our calculations. 

In summary, the interplay of Kondo and ferromagnetic 
(Stoner) interactions in large quantum dots leads to in- 
teresting ground state spin properties. The transitions to 
higher spin ground states due to the Stoner interaction 
are shifted non-monotonically in the presence of a Kondo 
spin. At weak Kondo coupling the Kondo spin acts just 
as an external field to assist ferromagnetic polarization. 
At strong coupling the system is described (again) by a 
universal Hamiltonian, but with a renormalized exchange 
constant for a reduced dot with one less level and one less 
electron. Ferromagnetic polarization can be either en- 
hanced or reduced in this limit, depending on how close 
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FIG. 3: (Color online) Spin transition curves for N = 12 at 
Zeeman field B and Kondo coupling Jk- The exchange cou- 
pling is (a) J s — and (b) J s /S = 0.825. Colored solid (black 
dotted) lines show the numerical (perturbative) solutions. 

the dot is to the bulk Stoner instability. The weak and 
strong coupling limits are described well by perturbation 
theory and our exact numerical solutions find a smooth 
behavior in the non-perturbative crossover region. 
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